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Abstract 

It is shown that partial incoherence, in the form of stochastic phase noise, 
of a Langmuir wave in an unmagnetized plasma gives rise to a Landau-type 
damping. Starting from the Zakharov equations, which describe the nonlinear 
interaction between Langmuir and ion-acoustic waves, a kinetic equation is 
derived for the plasmons by introducing the Wigner-Moyal transform of the 
complex Langmuir wave field. This equation is then used to analyze the sta- 
bility properties of small perturbations on a stationary solution consisting of 
a constant amplitude wave with stochastic phase noise. The concomitant dis- 
persion relation exhibits the phenomenon of Landau-like damping. However, 
this damping differs from the classical Landau damping in which a Langmuir 
wave, interacting with the plasma electrons, loses energy. In the present pro- 
cess, the damping is non-dissipative and is caused by the resonant interaction 
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between an instantaneously-produced disturbance, due to the parametric in- 
teractions, and a partially incoherent Langmuir wave, which can be considered 
as a quasi-particle composed of an ensemble of partially incoherent plasmons. 
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More than fifty years ago, Landau |l] discovered a damping of electron plasma waves 
whose linear dynamics is governed by the Vlasov-Poisson system of equations. This damping 
was caused by resonant wave-electron interactions. The Landau theory clearly shows that 
the decay rate of the wave energy due to the interaction with the resonant electrons of 
the system (satisfying p = u/k, where p denotes the electron velocity and u and k are 
the frequency and wave number of the electron plasma wave ) is proportional to the first 
derivative of the equilibrium distribution function po(p) of the electrons. Typically the shape 
of po(p) is such that dpo(p = u/k)/dp < 0, which implies that there are more particles with 
p < uj/k (which gain energy from the wave) than with p > u/k (which give energy to the 
wave). Consequently, the net result is that the wave is damped. However, finite amplitude 
Langmuir waves in plasmas can be created when some free energy sources, such as electron 
and laser beams, are available in the system. In the past, several authors [0- |J have 
considered the nonlinear coupling between high-frequency Langmuir and low-frequency ion- 
acoustic waves. Physically, this coupling occurs because a large amplitude plasma wave, 
interacting with small density ripples, produces a current which becomes the source for an 
envelope of Langmuir waves. The low-frequency ponderomotive force of the latter reinforces 
the density modulation in such a way as to produce a space-charge electric field, which 
varies on a time scale longer than the electron plasma period. Thus, within the slowly- 
varying envelope approximation, the dynamics of the Langmuir wavepacket is governed by a 
nonlinear Schrodinger equation (NLSE), while the non-resonant density ripples follow an ion- 
acoustic wave equation driven by a ponderomotive force created by the Langmuir wave). In 
his classic paper, Zakharov || demonstrated that these nonlinearly coupled equations admit 
a class of modulationally unstable stationary solutions and also exhibit the phenomenon of 
wave collapse. 

In this paper, we present a Landau damping of partially coherent Langmuir waves whose 
dynamics is governed by the Zakharov equations. In fact, we show that the latter can be 
reduced to a kinetic-like equation for a Wigner-Moyal-like transform for the amplitude of 
the electric field associated with the Langmuir wave. Based on this kinetic equation we 
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consider the dynamics of small perturbations on a partially coherent background solution 
consisting of a constant amplitude and a stochastically varying phase. The concomitant 
dispersion relation depends strongly on the spectral energy density of the Langmuir wave 
and a Landau-type damping of the perturbations is found due to the resonant interaction 
between the Langmuir quasi-particle and the instantaneously created ripples. 

As it is well known, the nonlinear propagation of a one-dimensional Langmuir wave- 
packet in an unmagnetized plasma with unperturbed density no, is governed by the Zakharov 
equations [|J 

BE 3vl d 2 E 5n „ , 
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where E is the slowly varying complex electric field amplitude, 5n is the density fluctuation, 
Vte is the electron thermal velocity, u pe is the electron plasma frequency, T e is the electron 
temperature, c s is the sound speed, and E is the unperturbed electric field amplitude. 
Moreover, x and t play the role of configurational space coordinate (longitudinal coordinate 
with respect to the wavepacket centre) and time, respectively. We will show that the nonlin- 
ear propagation of a Langmuir wavepacket, governed by the above Zakharov system, can be 
suitably described in phase space by a kinetic-like equation within a framework similar to 
the one of the Vlasov-Poisson system. A similar approach has recently been used to study 
the modulational instability of a large-amplitude electromagnetic wavepacket as well as the 
coherent instability analysis of intense charged-particle beam dynamics Within this 

framework, we show that a phenomenon similar to the classical Landau damping occurs 
for partially coherent plasma waves. In particular, this phenomenon tends to suppress the 
modulational instability. To this end, the above system of equations Q0.1] ) and (|0.2| ) can be 
cast in the following form 

a 2 d 2 ^ 
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where the angle brackets account for the statistical ensemble average dla Klimomtovich || 
due to the partial incoherence of the waves, s = ^/SXpe^pet (^De = v te /to pe is the electron 



Debye length), a = \/3Ao e , U = 5n/2n , /j, = ym e /3mj (m e and are the electron and 
the ion masses, respectively), and 

* = -7^— , (0-5) 
^8im T e 

= 'S^ . (0.6) 



Note that flCTq ) can be written as 

po = |^o| 2 = M 2 — , (0.7) 

where e = Eq/8tt and = raoT e are the electric energy density and the electron thermal 
energy density, respectively. According to the theory modelled by the Zakharov system 
( p.3p and (p.4[), e/e^ >> 1 (i.e, the ponderomotive effect is larger than the thermal pressure 
effect). Note also that Eq. ( p.4p implies that U is a functional of \^\ 2 , namely U = 
U [\^\ 2 ]. Consequently, once (|0.4j) is combined with (|0.3| ), the latter becomes a nonlinear 



Schrodinger equation (NLSE) with the nonlinear term [/[l^l 2 ]. This implies that, as in 
Quantum Mechanics, one can introduce the density matrix | as a sort of "two-points 
correlation function" and then the Wigner-Moyal transform which is solution of the 
von Neumann- Weyl equation In case of surface gravity waves, whose dynamics and 



instability is described in deep water by a nonlinear Schrodinger equation, such a kind of 
approach has been given in the pioneering works by Alber ||12|| , Crawford et al. [13[] and 



Janssen [f[3| |. More recently, a similar approach has been developed for the propagation of 
electromagnetic wavepackets in nonlinear media where an appropriate kinetic equation is 
able to show a random version of the modulational instability |B|.0,|IJJ. In particular, this 
approach can be also seen as an application of the Klimontovich statistical average method 



15| . The transition to the phase-space x—p (p = dx/ds is the conjugate variable of x) allows 
us to write a von Neumann equation || for w(x,p, s). In fact, although, we are dealing with 
a classical system, by following Alber [12], one can introduce the correlation function (whose 



corresponding meaning in Quantum Mechanics is nothing but the density matrix for mixed 

states) 

p(x,x',s) = (V*(x,s)y(x',s)) , (0.8) 

where the statistical ensemble average takes into account the incoherency of the waves. 
The technique was successfully introduced in the field of statistical quantum mechanics to 



describe the dynamics of a system in the classical space language [fLzfl - [14fl . Thus one can 
define the following Wigner-Weyl transform || 

w ^ s)= ^LSl p { x+ l x ~l s ) exp G?) dy ' (o - 9) 

which allows us to transit from the configuration space to the phase space. In fact, w(x,p, s) 
is governed by the following kinetic-like equation (von Neumann-like equation || ) 

dw dw _~ (-l) n fa\ 2n d 2n+1 Ud 2n+1 w 

+ P ~dx~ ~ ^ (2n + 1)! V 2 ) dx 2 ^ 1 dp 2 ^ 1 ' ^ ' ' 

On the other hand, since ( p.9|) implies that 

/CO 
w(x,p,s)dp , (0.11) 
-00 

the functional [/[l^l 2 ] can be expressed as a functional of w, namely U = U[w]. Conse- 
quently, (P-4|) can be written as 

(d 2 d 2 \ d 2 f r°° r°° \ 

a? - " a?J u = "a? (L m *" L m ° *) • (0 ' 12) 

where wq is the Wigner-Moyal transform of (po = |^o| 2 — S^°oo w o dp). The above scheme 
allows us to carry out a stability analysis of the stationary Langmuir wave by considering the 
linear dispersion relation of small amplitude perturbations. We start from the equilibrium 
state: w = w (p), U — Uq — 0, and perturb the system according to 
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w(x,p,s) =w (p) +Wi(x,p,s) , (0.13) 



U(x,s) = U + = C/i(x,s) , (0.14) 



where u>i(x,p, s) and C/i(x, s) are first-order quantities. Consequently, ( p.lOp and ( p.!2| ) can 
be linearized as follows 

dw! dw x ^ (-1)" (a^ cP^Hh (2n+1) 

da +P dx ~ £i(2n + l)\\2j dx^ W ° ' (U ' i5j 

/ Pp Pp \ fp / r oo \ 

^ )Ui = -i-A / , (0.16) 



y ds 2 dx 2 J dx 2 

where w^ n+1 ^ = d 2n+1 po/dp 2n+1 . By introducing the Fourier transform of Ui(x,s) and 
u>i(ac,p, s), i.e. 

/oo /.oo . 

d/c / dcu Ui(k,uj) exp (ikx — iujs) , (0-17) 
-oo J —oo 



oo 



Wi(x,p,s)= / dk du wi(k,p, uj) exp {ikx — ius) , (0.18) 



and substituting ( p.!7|) and ( |0.18| ) into Eq.s ( p.!5| ) and ( |0.16| ), we readily obtain the following 
dispersion relation 

ji _ 2 k 2 = k 2 f°° w (p + ak/2) - w (p-ak/2) dp 

J-oo ak p — uj/k 

We consider a Lorenzian spectrum of the form: 



MP) = (0-20) 
7i A z + p A 



where A is the width of the spectrum. Dispersion relation ( p.!9| ) becomes 



(u 2 - fi 2 k 2 ) (u 2 - a 2 k 4 /A + 2ikAu - k 2 A 2 ) = p k 4 , (0.21) 

As the simplest case, we take the limit of A — > 0. This case corresponds to a coherent 
monochromatic wave-train and the dispersion relation reduces to 



(uj 2 -fi 2 k 2 ) (uj 2 -a 2 k A /A) = Po k A , (0.22) 

which predicts the well known result of a purely growing instability M . 

In the following, we concentrate on the limit of small ak << 1 and u/k » 1. We 
will show that during its nonlinear evolution the Langmuir wavepacket can exhibit a phe- 
nomenon, similar to the Landau damping [PJ, between parametrically-driven non-resonant 
density ripples and Langmuir quasi-particles. It is easy to verify that the condition ak « 1 
(i.e., \oek « 1) implies that 

w (p + ak/2) — w (p — ak/2) 



ak 

Consequently, (|0.19| ) becomes 



dwo/dp = w . (0.23) 



c 2 = k 2 I dp , (0.24) 

J p — u/k 

where, using the assumption u/k » 1, we can neglect the term involving /i 2 because 
/i << 1. 

Note that ( p. 24] ) is formally identical to the linear dispersion relation that holds for warm 



plasma waves |T[ or for non-monochromatic charged-particle bunches in circular accelerating 
machines for an arbitrary complex coupling impedance fl7| . Consequently, the dispersion 



relation (|0.24j) allows us to predict a sort of Landau damping Jy], as described in plasma 
physics for linear plasma waves as well as in charged-particle beam physics. A similar effect, 
which has been called quantum-like Landau damping (QLLD), has recently been predicted 
for the nonlinear propagation of a large-amplitude electromagnetic waves in nonlinear media, 
such as optical fibers and plasmas as well as for the quantum-like description of charged- 
particle beams in an accelerating machine [|6|,|7|,|T5f . In order to investigate this phenomenon 
also for large-amplitude Langmuir waves, we apply the standard procedure [H| used for the 



Landau damping of a linear plasma wave, which involves integrating the above dispersion 
relation in the complex plane, evaluating both the Cauchy principle value (PV) and the 
semi-residue in p = u/k. The result is 
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u 2 = k 2 I W ° - dp + mk 2 w (f3) , (0.25) 
Jpv p — uj/k 

For a symmetric sufficiently smooth stationary background distribution, Wq(p), we have 

/ 3A 2 ^ 



,j 2 ( 1 + — j k 2 + i7ck 2 w' {(3) , (0.26) 



where A = ((p 2 )) 1 ^ 2 denotes the r.m.s. width of w (p), and (3 = p\j 4 (1 + 3A 2 /4pg). Note 
that, according to the above hypothesis, we have (3 » 1 and A 2 /p << 1- This implies 
that p = |\l>o| 2 >> 1, and from (|077| ) it follows that e/e T » 1//J 2 . Eq. (|0.26|) clearly 



shows that the damping rate is proportional to the derivative of the Wigner distribution u;o- 
This is formally similar to the expression for the Landau damping of a linear plasma wave 
in a warm unmagnetized plasma. In fact, writing u in the complex form u = ujr + i&j, 
substituting it in (|0.25|) and then separating the real and imaginary parts, we obtain 



o*k - = k 2 I -^jr dp , (0.27) 
Jpv p — oo/ k 



and 



ui = 1—™'M ■ (0-28) 



The Landau damping rate of the Langmuir wave in the appropriate unities is 7 
^ujpeXceUJi. For example, for a Gaussian wavepacket spectrum, i.e. wo(p) 
(p /V27rA 2 ) exp (— p 2 /2A 2 ) } one obtains 



3tt \^ 2 p 5 ' 2 f 6 y/ 4 



7 = - ( — ) -^3- y—j oj pe X De k exp 



/i / r \ 1/2' 



2 A 2 Ve r 



(0.29) 



where higher-order terms have been neglected. The present damping mechanism differs from 
the standard Landau damping in that here Wq(p) does not represent the equilibrium velocity 
distribution of the plasma electrons, but it can be considered as the "kinetic" distribution 
of all Fourier components of the partially incoherent large-amplitude Langmuir wave in the 
warm plasma. In terms of plasmons, we realize that Wo(p) represents the distribution of the 
partially incoherent plasmons that are distributed in p-space (i.e. in fc-space) with a finite 
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"temperature". Consequently, the QLLD described here is due to the partial incoherence of 
the wave which corresponds to a finite-width Wigner distribution of the plasmons. 

In this paper, we have carried out an analysis for a Landau-type damping in plasmas 
containing partially incoherent Langmuir waves. Starting from the Zakharov system of equa- 
tions and employing the Wigner-Moyal quasi-distribution function for the complex Langmuir 
wave electric field amplitude, it has been shown that the Zakharov system of equations can 
be converted into a pair of coupled evolution equations in phase space. For a monochromatic 
and coherent wave (wo(p) = Po8(p)), we have obtained a dispersion relation which exhibits 
a purely growing instability |4]]. Furthermore, for a broad-band Langmuir wave spectrum, 
corresponding to a partially incoherent constant amplitude wave, we have investigated the 
properties of small perturbations for ak << 1 (i.e., \oek << 1), and very large dimension- 
less phase velocity (i.e., uj/k » 1). Our dispersion relation for this case exhibits a sort 
of weak Landau damping, very similar to the one predicted for the Vlasov-Poisson system 
for linear plasma waves. The physical origin of this phenomenon is attributed to the "non- 
monochromatic" behaviour of the Wigner spectrum of the Langmuir wave, which can be 
considered as an "ensemble of incoherent plasmons". Actually, a broad-band spectrum of 
incoherent plasmons, forming the "wave packet", interact individually with the self-excited 
density ripples. Similar to the standard Landau damping, where the electrons interact indi- 
vidually with a linear plasma wave and statistically produce a net transfer of energy from 
the wave to the particles, the plasmons interact individually with the ripples and produce 
a change of energy which is more significant around the resonance, which is determined by 
the condition p = u/k. Finally, we observe that since the system can be modulationally 
unstable this phenomenon acts in competition with the modulational instability. 
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